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""^3 Abstract 

We present the first purely semiclassical calculation of the resonance spec- 
■ trum in the Diamagnetic Kepler problem (DKP), a hydrogen atom in a con- 

stant magnetic field with L z = 0. The classical system is unbound and com- 
,-£h pletely chaotic for a scaled energy e ~ EB~ 2 / 3 larger than a critical value 

e c > 0. The quantum mechanical resonances can in semiclassical approxima- 
tion be expressed as the zeros of the semiclassical zeta function, a product 
over all the periodic orbits of the underlying classical dynamics. Intermit- 
tency originating from the asymptotically separable limit of the potential 
at large electron-nucleus distance causes divergences in the periodic orbit 
formula. Using a regularisation technique introduced in [G. Tanner and D. 
Wintgen, Phys. Rev. Lett. 75, 2928 (1995)] together with a modified cycle 
expansion, we calculate semiclassical resonances, both position and width, 
which are in good agreement with quantum mechanical results obtained by 
the method of complex rotation. The method also provides good estimates 
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for the bound state spectrum obtained here from the classical dynamics of a 
scattering system. A quasi Einstein-Brillouin-Keller (QEBK) quantisation 
is derived that allows for a description of the spectrum in terms of approxi- 
mate quantum numbers and yields the correct asymptotic behaviour of the 
Rydberg-like series converging towards the different Landau thresholds. 

PACS numbers: 0545, 0365.Sq, 31.15.Gy, 32.30 



The hydrogen atom in a uniform magnetic field has become one of the most im- 
portant examples for studying the correspondence between quantum mechanics and 
classical chaos. The Hamiltonian is known to a high accuracy and furthermore, the 
system is experimentally accessible in the laboratory. The angular momentum in 
the direction of the magnetic field is conserved, which reduces the classical system 
to a problem with two degrees of freedom. The classical flow in phase space covers 
a wide range of Hamiltonian dynamics reaching from bound, nearly integrable be- 
haviour to completely chaotic and unbound motion by varying one parameter, the 
scaled energy e. 

In the 1980's, the system served as a catalyst for quantum chaos. Modulations 
in the absorption spectra of highly excited hydrogen atoms in a magnetic field, the 
so-called quasi Landau levels [1, 2], could be assigned to classical trajectories [3, 4] 
and could be understood in the framework of the semiclassical periodic orbit theory 
[5, 6, 7]. This strikingly simple answer to a long outstanding question led to a 
completely new viewpoint in analysing quantum spectra in general. We refer the 
reader to the review articles and article collections [8, 9, 10, 11, 12]. 

So far, the main interest in the literature laid on understanding quantum man- 
ifestations of the system such as the quantum spectra or the wave functions, in 
terms of the underlying classical mechanics. Effective numerical techniques to 
solve Schrodinger's equation directly (neglecting relativistic and finite mass ef- 
fects) have been developed in the recent years, both for the bound state spectra 
[13, 14, 15, 16, 17] and resonances in the continuum [18, 19]. Periodic or closed 
orbits appear in the Fourier transformed energy - or photo ionisation - spectra ob- 
tained from experiments or from full quantum calculations. Classical orbits can in 
some cases be identified as "scars" in quantum eigenfunctions [20] and as prominent 
peaks in the corresponding Husimi distributions [21]. 

Less effort has been undertaken in the opposite direction, i.e. in quantising hy- 
drogen in a magnetic field directly in terms of classical entities only. The density of 
quantum states [22], photo absorption spectra [23, 24] or the evolution of Rydberg 
wave packets [25] can formally be written as sum over periodic or closed orbits. 
A brute force application of the orbit sums using the shortest classical orbits only 
yields the coarse grained structure of the spectrum [26] or the time evolution of 
Rydberg wave packets on short time scales [25], respectively. To resolve the fine 
structure of the spectrum and individual highly excited quantum states, the infor- 
mation of the classical dynamics on long time scales is needed. The classical orbit 
sums, however, diverge due to the exponential growth of the number of orbits in 
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chaotic systems and appropriate resummation techniques have to be developed to 
overcome this problem. 

Starting point is here the spectral determinant D(E) = Y\ n (E — E n ), which can 
in semiclassical approximation be written as product over all periodic orbits of the 
system, the so-called semiclassical zeta function [27]. An analytic continuation of 
the product formula may be given by expanding the product and regrouping terms 
with the help of a symbolic description of the flow in form of a cycle expansion 
[28, 29]. This method has been shown to work for strictly hyperbolic classical 
dynamics, where all periodic orbits are unstable and the Lyapunov exponents of 
the orbits, i.e. the logarithm of the largest eigenvalue of the Monodromy matrix 
divided by the period, are strictly bound away from zero. 

The Diamagnetic Kepler problem (DKP) corresponding to conserved angular 
momentum L z = 0, (with the magnetic field in z direction), is known to be com- 
pletely chaotic in the sense that all periodic orbits are unstable for values of the 
scaled energy e > e c = 0.328782 . . . [30]. The classical motion is unbound in this 
regime and the flow in phase space can be described in form of a complete ternary 
symbolic dynamics [31, 30]. The DKP is, however, not hyperbolic even for energies 
e > e c . The potential becomes separable in the limit where the classical electron is 
far from the nucleus, causing regular, but unstable, dynamics in this phase space 
region. As a consequence, the Lyapunov exponents of periodic orbits extending 
far into the regular region of phase space tend to zero. This behaviour, known as 
intermittency in the chaos-literature, is typically for the neighbourhood of marginal 
stable orbits or stable islands and is thus generic in Hamiltonian systems. 

The regular limit of the classical motion causes divergences in the semiclassical 
periodic orbit expressions. This kind of divergence is different from the well known 
convergence problem of semiclassical expressions due to the exponential growth of 
the number of periodic orbits. Divergences introduced through intermittency can 
be regularised by first summing over periodic orbit contributions in the regular 
classical regime alone. An analytic continuation of each of the various sums can be 
given explicitly [32]. This procedure is equivalent to a cycle expansion of the zeta 
function in terms of an infinite alphabet [29]. We can finally present the first purely 
semiclassical quantisation of the resonance spectrum of DKP at fixed e > e c . The 
regular limit of the dynamics can be identify to cause the Rydberg series structure 
at the various Landau thresholds. The high excited Rydberg resonances can be well 
described by a modified Einstein-Brillouin-Keller quantisation. 

The article is written self-consistently in such a way that it contains all the 
information needed to understand the classical, the quantum and the semiclassical 
aspects of the problem. The classical system together with the symbolic dynamics 
is introduced in section 1. Asymptotic expressions for the actions and stability 
exponents in the separable limit of the potential are derived here. In section 2, 
we explain the method of complex rotation, which is an efficient tool to treat the 
full quantum problem in cases of resonances in the continuum. In the last section, 
we present our new semiclassical quantisation, and results are compared with full 
quantum calculations. 
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1 Classical dynamics 



The non relativistic classical Hamiltonian for the hydrogen atom in a uniform mag- 
netic field with field strength B along the z-axis is given as 

H = 2^ " 7 + \ m ^ 2{ - x2 + ^ + w 

The 2 - component of the angular momentum L z is conserved and we will restrict 
ourself to the problem L z = in the following. We work in the infinite nucleus 
mass approximation and m e denotes the mass of the electron. The frequency uj = 
eB/2m e c is half the cyclotron frequency. The Hamiltonian (1) in atomic units and 
for L z = has the form 



H = E = Ei + i- 1 + l 7 y (2) 

2 2 vV 2 + 2 2 8 r F V 7 

with the magnetic field strength 7 = B/B written in units of B = mle 3 c/h = 
2.35 • 10 5 T. The radial distance from the z-axis is given by the coordinate p. The 
system is bounded for E < 0, while for E > almost all trajectories escape to 
z = ±00 with nonzero kinetic energy. 

Introducing the scaling transformation 

r = 7 - 2 / 3 ?, p = 7 1 /3p ) (3) 
yields the new Hamiltonian 

iJ = £ = 7 - 2 /3 W = f + f__L_ + ^, (4) 

which is independent of 7. The classical dynamics is now controlled by one param- 
eter only, the scaled energy e = E^~ 2 I Z . The classical action along a trajectory 
scales with the magnetic field like 



J pdq 
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-V 3 S (5) 



with e fixed. 

The Hamiltonian (4) is singular at f = yielding singular equations of motions 
at this point. For solving the equations of motion numerically, it is more convenient 
to switch to semi-parabolic coordinates and momenta [8] 

v 2 = f — z, /j, 2 = f + z, Pv = dv/dr, = d/i/dr, (6) 
with a new rescaled time r 

dr = dt/2f = {v 2 + ii 2 y l dt. (7) 
The Hamiltonian (4) in semi-parabolic coordinates is 

h = \pl + \p\ - <v 2 + // 2 ) + ^ VV + /i 2 ) = 2, (8) 

and the scaled energy e enters as a parameter here. 

The structure of the dynamics of the Hamiltonian (8) depends on the value of 
the scaled energy e alone. The E-^ parameter plane can be partitioned into 5 
distinct regions [8, 30], see also Fig. 1: 
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Figure 1: The 5 main dynamical regions in the E-^y parameter space reaching from 
bound and almost integrable dynamics (I) to unbound completely chaotic dynamics 
(V). See also the explanation in the text. 

(I) e < —0.5: bounded almost integrable motion; 

(II) —0.5 < e < —0.13: bounded motion with mixed chaotic and regular motion; 

(III) —0.13 < e < 0.0: the last large stable island disappear; the dynamics is mostly 
chaotic; 

(IV) 0.0 < e < 0.328782...: unbounded mostly chaotic motion; the symbolic 
dynamic is not complete; 

(V) 0.328782 . . . = e c < e: unbounded chaotic motion with a complete symbolic 
description. 

In the following, we will always work in the region e > e c = 0.32878 . . ., where 
a simple and complete symbolic description can be assigned to flow in phase space 
[30]. Furthermore, all periodic orbits are unstable in this parameter regime. For e 
values below the critical energy e c , the symbolic dynamics become incomplete and 
stable islands appear. Stable classical motion implies additional complications in 
the semiclassical description in section 3, which will not be included in the work 
presented here. 

1.1 Symbolic dynamics and periodic orbits 

We can motivate the structures of the classical motion and the underlying sym- 
bolic description by first studying the dynamics of a symmetric four disk scattering 
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billiard [33]. The similarity between the DKP and the four disk system becomes 
obvious, when looking at the shape of the potential in the regularized fi, v - coor- 
dinates, see Fig. 2 (b). In the four disk billiard, a trajectory may after a bounce 
with one disk either hit one of the other three disk or escape from the system. A 
trajectory bouncing two times after leaving the first disk may bounce in 9 different 
ways. We find that there are 3" possible sequences of disks for n bounces. This set 
of non-escaping starting points for n — >■ oo is called a Cantor set. We label each 
disk with a number s t € {1, 2, 3, 4} and a string Si ■ ■ ■ s n assigned to a trajectory 
represents one of the 3" possible sequences. A particle cannot bounce twice off the 
same disk, which is the only restriction in the 4-letter symbolic dynamics as long 
as the distance between the disks is sufficiently large compared to the disk radius. 
The orbits remaining in the system both forward and backward infinitely long in 
time form a two dimensional Cantor set and are uniquely described by bi-infinite 
symbol strings • • • s_ 2 s_iS SiS2 ' ' '• 

Since the Cantor set is ternary we can reduce the number of symbols to three 
with no restrictions on the sequence of symbols. The reduction can be done in 
several possible ways [31, 33]. In Appendix A, we construct a ternary alphabet 
closely related to the C± v - symmetry of our system, which will be used throughout 
the paper. A trajectory is now characterized by a symbol string • • • g-igogi ■ ■ ■ with 
gi £ {0, 1, 2} defined in Appendix A. 

In Appendix B, we introduce a so called well ordered symbolic dynamics • • • w_iWqWi ■ • • 
with Wi € {0, 1, 2}, which is in particular useful for finding periodic orbits numeri- 
cally [34]. 

The structure of the Cantor set is similar for the DKP and for the four-disk 
scattering system. We can associate a 4-letter symbol string to every trajectory 
also in the DKP, and the symbols can be understood as a smooth "bounce" with one 
of the steep hills in the potential (8). At this bounces, a trajectory has a caustic, 
and we can determine the symbolic dynamics either by determining the number of 
caustics [31] or by making a suitable partition curve in a Poincare plane [35]. A 
detailed discussion of this method will be given elsewhere [36]. 

The symbolic dynamics is complete, i.e. each possible symbol string • • • g^ig^gi • • • 
corresponds to a non-escaping orbit in the symmetric four disk system, if and only 
if the gap between two disks are larger than 0.205 times the disk radius[37]. If 
the disks come closer to each other certain paths between the disks are forbidden. 
These paths lie now in the "shadow" of a disk in between. The corresponding sym- 
bol strings cannot be related to a physical orbit and the symbolic dynamics is said 
to be pruned [38, 39, 37]. 

We find an analogous situation in the DKP. Here, each possible symbol string 
corresponds to an orbit if the scaled energy is larger than e c = 0.32878 . . . [30]. For 
energies below the critical value, the hills of the potential will shadow some of the 
trajectories which exist for larger e and the symbolic dynamics becomes pruned. 
The value of e thus plays in the DKP the same role as the distance between the 
disks in the 4-disk system. Note that pruning in the DKP is introduced through 
bifurcations and stable islands appear and disappear with varying e < e c . In the disk 
system, however, periodic orbits are always unstable and disappear immediately at 
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Figure 2: Periodic orbits in the (a) (z, p) plane and (b) (u, p) plane. The full line 
denotes the boundary of the potential for e = 0.5. 

a bifurcation point. 

The periodic orbits are a subset of all non-escaping orbits existing in the system, 
(i.e. of the repellor). A periodic orbit is described by an infinite repetition of a finite 
symbol string G = gig 2 ■ ■ ■ g n - In the following, we focus mainly on periodic orbits, 
which we characterize by the shortest non repeating symbol string directly. 

Four short periodic orbits of the DKP are drawn in Fig. 2, both in the coordinate 
space (z,p), and in the semi parabolic coordinates (i>, p). Note that the two orbits 
labelled G=2 and G=20 in Fig. 2 (a) collide with the nucleus. In Fig. 2 (b) we 
have indicated the labelling of the four hills used to construct the 4-letter symbolic 
dynamics. A translation from the 4 - letter to the 3 - letter alphabet can now be 
read off directly with the help of Appendix A. The orbits labelled 1, 2, 10, and 
20 in Fig. 2 (a) are given the symbolic descriptions S = 1234, 13, 1232, and 1324, 
respectively, in the 4 - letter code, Fig. 2 (b). 

1.2 Asymptotic in the classical Landau channel 

In this section, we will study the dynamics of the Hamiltonian (2) for positive energy 
E in the limit \z\ — > oo. The coupling between the z and p coordinates vanishes in 
that limit and the Hamiltonian (2) has the asymptotic form 

H = \ P 2 P + \vl + y P 2 - ± + Otf/z*) = E > 0, (9) 

with lo = 7/2. The motion separates for large \z\ values in a pure Coulomb part 
along the z - axis and an harmonic oscillator perpendicular to the magnetic field. 
The energies 

E h = \p 2 p + Y p2 and Ec = l p2z + U (10) 
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become adiabatic constants of motion for \z\/p —> oo. The regular region far from 
the nucleus is always coupled to the chaotic motion near the core due to the at- 
tractive Coulomb force in z - direction. The classical motion of the electron in 
the DKP thus alternates between strong chaotic motion for p m \z\ and regular 
time intervals out in the channel \z\ p. Almost all trajectories escape finally 
to infinity with nonzero momentum in z - direction. Orbits which go further and 
further out in the regular channel but return to the origin, approach a marginally 
stable periodic orbit \z\ = oo, p z = 0. In the symbolic description explained in the 
appendix A, this marginally stable orbit is labelled 0. Escaping trajectories which 
continue travelling along the z - axis can be separated from returning once by a 
plane, that is part of the stable manifold of the marginal stable orbit. This plane 
is in the asymptotic limit given by the condition E c = p 2 z /2 — l/\z\ = 0. Returning 
trajectories can be associated with bound motion in the z - direction and E c < 0. 

Periodic orbits that go far out in the regular region \z\ ^> p pick up regular 
contributions to the actions and stability exponents, which can be given in analytic 
form to leading order. In the ternary symbolic description the channel orbits are 
characterized by a long string of consecutive symbols '0' in the symbol code. Of 
particular interest are periodic orbits with a symbol code of the form GO 11 , i.e. 
periodic orbits with a common head string G and a tail of n symbols '0'. We will 
call these orbits a periodic orbit family. The orbits in a family have approximately 
the same behavior before entering and after leaving the regular region. They differ 
only in the '0' - tail string, i.e. in the number n of half - oscillations perpendicular 
to the z - axis in the regular region. In Fig. 3, some members of the periodic orbit 
families with head string G = 1 are shown explicitly. 

The leading n - dependent behavior in the actions and stability exponents is 
universal, i.e. it does not depend on the past or future of the periodic orbits before 
entering and after leaving the regular phase. Non universal contribution approaches 
a constant for orbits with the same head string G, but increasing number of oscil- 
lations in the regular regime. 

To obtain the leading contributions to the actions, we start with the Hamiltonian 
(9) written in action-angle variables (Jh, (fh, J a fc), 



The integration is taken over a single revolution of a trajectory in the potential of 
the one dimensional Hamiltonians in (10). (Note, that p > in (12)). The angle 
variables propagate linearly in time with frequencies 



H S ep — %^Jh — 

and omitting the coupling terms. The action variables are given as 



(11) 




(12) 



(13) 



u h = 




2co; oo c = 




(14) 
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Figure 3: Three members of the periodic orbit family 10" with (a) n = 25, (b) 
n = 50, and (c) n = 100 in the fundamental domain fi > v > of the regularized 
(i, v coordinates. 
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The surface of constant J^, J c form a two dimensional torus in the 4-dimensional 
phase space. The total action of a trajectory on a torus after the time T c = 2tt /lu c , 
i.e. after one oscillation in z - direction, is 

J(a,E) = J c (a,E) + aJ h (a,E), (15) 

where the winding number a = LOh/ui c is given as 

a = 2u.Jl (16) 

Inserting (12), (13) and (16) in (15), we can write the total action as function of 
the energy E and en, 

3 ^ E) = tA e+ ^) \ (17) 

In terms of the scaled energy e = E^ 2 ^ and 7 = 2co, we obtain 

J(«,e)=« 7 - 1 /3^ + ^-2/3^ (18) 

The symbol '0' corresponds to one oscillation perpendicular to the field axis in the 
p coordinates or half an oscillation in the original (x, y) coordinates. A trajectory 
with n consecutive zeros in the symbol code has an approximate winding number 

a(n) m n + a for n ^> 1. (19) 

Contributions due to the non - integrable part of the dynamics are collected here 
in the shift term a , which approaches a constant in the limit n — >■ 00 for periodic 
orbits in the same family. Inserting (19) in (18), the asymptotic n dependent 
behavior of the (un-scaled) action of a channel trajectory can be written as 



S{n) = 27T7" 



-1/3 



n 



(e + ^ 2/3 ) + So + 5 2/3 n" 2 / 3 + ^(n- 1 )] . (20) 



The coefficients S , S2/3, ■ ■ ■ depend on e and explicitly on the head string G or on 
the past and future of the orbit before entering and after leaving the regular phase. 
Fig. 4 (a) shows results for the periodic orbit families 10" and 20 n up to n = 100. 
The constant family dependent term in (20) is So = —0.6483 for the 10" family and 
S = -0.1440 for the 20" family, (see also Table 1). 

The asymptotic behavior of the stability exponents for periodic orbits with long 
sequences of symbols '0' can be extracted from the trace of the stability- or Jacobi- 
matrix of the classical flow, 

d(qo,Po) 

This matrix describes the evolution of the classical dynamics in the neighborhood 
of a given trajectory (q(t),p(t)) in linear approximation. 
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Figure 4: The actions (a) and stability exponents (b) of the periodic orbit families 
10™ (O) and 20™ (+) plotted with respect to the number n of zeros in the symbol 
code (for e = 0.5); the full lines correspond to the asymptotic expressions S n ~ 
n(e + 3/2n~ 2//3 ) and X n = 5/31og(n), respectively. 



The asymptotic behavior of TrM for a large number n of symbols '0' is derived 
in Appendix C and we state only the final result here, 

TrM = /cin 5/3 + k 2 n + 0(n 2/3 ). (22) 

The terms ki,k 2 depend on e and on the past and future of the periodic orbit 
before entering the regular region. The largest eigenvalue of the stability matrix 
grows approximately like the trace of M. The stability exponent of a periodic 
orbit, i.e. the logarithm of the largest eigenvalue of the stability matrix, scales thus 
asymptotically like 

X(n) = | logn + l + l 2/3 n-s + C^n" 1 ). (23) 

The logarithmic behavior and the factor 5/3 is universal, i.e. independent of the 
scaled energy e and of the head string G. In Fig. 4 (b), we compare periodic orbit 
data of the families 10™ and 20™ with the leading term 5/3 logn. The curves deviate 
not more than a constant in the large n limit. This constant is again family depen- 
dent and we obtain l = 1.696 for the 10™ - family and l = 3.675 for the 20™-family 
(see Table 1). This is the slow logarithmic increase in the stability exponents, which 
reflects the regularity of the dynamic in the separable limit of the potential. The 
DKP shows the typical behavior of intermittency, i.e. the the Lyapunov exponent 
A = X(n)/T(n) ->■ for n ->■ oo. 

The discussion above is strictly valid only for scaled energies e larger than the 
critical value e c . The symbolic dynamics becomes pruned for e < e c and long strings 
of '0's tend to be forbidden below the critical value. The actual pruning rules depend 



11 



on the past and future of the '0' code strings. The maximal possible length of the 
'0' string in a periodic orbit family, which becomes pruned, decreases with e. The 
10" family is an example showing this behavior. There are also families, as e.g. 
the 20" family, where the limit n — >■ oo exists for all e > 0. A detailed analysis of 
these pruning rules will be presented elsewhere [36]. Note, however, that periodic 
orbit families with an infinitely long '0' tail allowed will have the same asymptotic 
behavior as discussed above. In the semiclassical analysis in section 3, we choose 
e = 0.5 to avoid problems due to the incomplete symbolic dynamics. 



2 The Quantum Spectrum 

The starting point of our quantum calculations is the Hamiltonian of a hydrogen 
atom in a uniform magnetic field, where we neglect relativistic corrections and 
effects of the finite nuclear mass recently discussed in [40, 41]. These effects occur on 
extremely long time scales, and correspondingly narrow energy spacings, which are 
accessible neither to present laboratory experiments nor to our numerical quantum 
calculations, and are also ignored in the classical and semiclassical calculations. The 
quantum Hamiltonian reads in atomic units and setting 7 = B/B 

H = \p 2 -1 + \^ + ^V, (24) 

in analogy with its classical counterpart (2). Exact quantum numbers of the Hamil- 
tonian (24) are the parity with respect to the (z = 0) - plane, n z , and the magnetic 
quantum number, m, but the system remains non separable in (p, z) coordinates. 
Eigenvalues and eigenstates of Schrodinger's equation can be calculated in each rrf 
subspace separately. 

Introducing dilated semiparabolic coordinates similar to (6) 

fj. = ^Vr + z ; v = ^Vr - z, (25) 
where b is a free length scale parameter, Schrodinger's equation is transformed to 

+ Aj, + 2£&V + » 2 ) " ~/l 2 A V + v 2 ) + 46 2 ] * = (26) 

where A^ and A„ are the Laplace operators for the radial coordinate in two dimen- 
sional space (p, <p) and (u, tp) with constant angular momentum L = m, i.e. 



19 d in 

A p = - — (p— ) 3-, (p = porz/), 

p op op p l 



The paramagnetic term \^L Z at constant magnetic field strength results in a con- 
stant energy shift = mj/2, which is already omitted in (26). 

The Schrodinger equation (26) depends on the two physical parameters E and 7. 
Eigenvalues can now be calculated along various cuts in this plane, e.g. at constant 
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magnetic field strength 7, along lines £7 = const or along lines of constant scaled 
energy e = E^~ 2 I Z . For example at constant magnetic field strength we introduce 

A = -(1 + 2E6 4 ), 

and obtain a generalized eigenvalue equation for the parameter A 

= A (/i 2 + v 2 ) * (27) 

with 

H = (A,-n 2 ) + (A u -v 2 ). (28) 

Due to scaling laws of the classical Hamiltonian (2), the structure of the clas- 
sical dynamics depends only on the scaled energy e = E^ 2 ^ . For comparison 
with semiclassical theories, it is therefore most convenient to evaluate the quantum 
spectra along lines e constant. Choosing 

s = 7~ 1//3 ; b = j3s = const 

with j5 a new free parameter, (26) can be transformed into a generalized eigenvalue 
expression for s 



* = s- 2 (Ap + A,) * . (29) 



Both Eqs. (27) and (29) can now be represented in matrix form as a generalized 
eigenvalue problem with sparse symmetric matrices in terms of the complete set of 
basis functions given by products of the eigenstates \N p m> of the two-dimensional 
harmonic oscillator (28), viz. 

\N li N v m> = \N^m> <g> \N u m> ; 7^,7^ = 0,1,2,... . (30) 

Eigenvalues A = — (1 + 2b 4 E) or s = 7 -1 / 3 and the corresponding eigenstates are 
obtained by numerical diagonalization of matrices. 

Because the basis states (30) are square integrable £ 2 functions scattering wave 
functions of open systems can not be expanded in this basis and the method is 
usually restricted to the calculation of bound states at energies below the ionisation 
threshold, i.e. 

E < E IP = ^{\m\ + 1) 

or for constant scaled energy e > 



7 -V3 < 



\m\ + 1 



2e 

Note, that the spectrum of hydrogen in a magnetic field exhibits bound states at 
energies even above the classical threshold E = 0, i.e. when the classical motion is 
already unbound in z direction. The continuum wave functions are (for \z\ — >■ 00) 
still quantized in the p - coordinate and the quantum threshold is thus shifted by 
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the zero point energy of the harmonic oscillator in the classical Landau channel. To 
account for unbound resonances in the continuum we adopted the complex-rotation 
method [42, 43], which is based on the replacement 

r — >re ie 

in the Hamiltonian and in the wave functions. By this transformation, hidden reso- 
nances of the Hamiltonian in the continuum, associated with complex eigenvalues, 
are exposed, while the resonance wave functions can still be described by the C% 
integrable basis functions (30), but with complex arguments. In our approach, the 
complex rotation by the angle 9 is reflected in the replacement of b and j3 in Eqs. 
(27) and (29) with the complex dilatation parameter 

b=\b\e i i and /5=|/5|e^, 

respectively. Representing Eqs. (27) and (29) with complex rotated coordinates in 
matrix form, this results in a generalized eigenvalue problem for complex symmet- 
ric, non-Hermitian matrices. Numerically, the eigenvalue equations were solved by 
extending the Lanczos algorithm [44] to complex matrices. 

If the scale parameter b or j5 and the complex rotation angle 9 are appropriately 
chosen and a sufficiently high number of basis functions (30) is considered in the 
numerical diagonalization isolated bound states along the real well as reso- 

nances in the lower complex half plane can be obtained to high numerical precision. 
The method fails close to the Landau thresholds 

S w=7 (jV + M±i) (31 ) 

or, for constant scaled energy e > 

SN = l(^ N+ H±3 j ,32) 

with N = 0, 1, 2, . . ., the Landau quantum number. To each Landau channel belongs 
an infinite number of Rydberg like bound states (N = 0) or resonances (N > 0) 
converging to the threshold. Obviously such infinite series of states can not be 
obtained from a calculation with a finite basis size. The numerical convergence of 
resonances close to the thresholds is extremely slow if all basis functions up to a 
given oscillator quantum number, i.e. 

Nn + N v + \m\ < iV max 

are considered in the calculation because the total basis size increases quadratically 
with growing N max . To improve the convergence we chose a different method to 
truncate the basis. In coordinate space states |A^ jU A^ iy m> with m N u are localized 
close to the (z = 0)-plane. In this direction the classical motion is bound and there 
is no need to complete the basis with states of this type. On the other hand states 
with <C N„ or N v <C are oriented along the z-axis, i.e. where the potential is 
open. These states are important for the representation of Rydberg like resonances 
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reaching far into the Landau channels. We have now optimized our basis by taking 
into account only those states which cover the physically important region. With 
this choice the numerical efficiency in the calculation of Rydberg like resonances 
could be considerably improved. 

Results for a quantization along constant e-lines will be presented in the next 
section in comparison with the semiclassical calculations. 

3 Semiclassical quantization 

In the following, we will present a method for a semiclassical quantization of hydro- 
gen in a constant magnetic field for L z = 0, i.e. with azimuthal quantum number 
m = 0. The quantization will be performed with respect to the magnetic field 
variable s = 7 -1 / 3 introduced in the last section along constant scaled energy lines 
e = £ 7 - 2 / 3 > e c . 

Starting point of our semiclassical calculation is the zeroth order Gutzwiller- 
Voros zeta function [27], which is the leading term in a semiclassical approximation 
of the spectral determinant det(E — H). The zeta function can be derived from 
the periodic orbit trace formula [22, 45] and has for the scale invariant classical 
Hamiltonian (2) with fixed scaled energy e the form 

CHs) = n (1 - t p (s)) , with t p (s) = e^p-iM-Tf. (33) 
v 

The product runs over single repeats of all periodic orbits of the system. The scaled 
action S introduced in (5) is taken here along a periodic orbit. The (complex) vari- 
able s is the scaling factor 7 -1 / 3 in units of Ti. The winding number a corresponds 
to the number of turns of the stable or unstable invariant manifold around the pe- 
riodic orbit [46]. The linear stability exponent A is defined as the logarithm of the 
largest eigenvalue of the stability matrix M, Eq. (21). (We restrict ourself here to 
the purely hyperbolic case, i.e. we presume that all periodic orbits are unstable and 
e > e c ). 

The semiclassical approximation to the quantum eigenvalues are given by the 
zeros of the zeta function. The magnetic field dependence enters the semiclassical 
zeta function for fixed e only through the scaling factor s = 7 ~ 1//3 . The winding 
number and the stability exponent are independent of 7 due to the scale invariance 
of the classical Hamiltonian. A quantization with respect to the magnetic field 
variable s uses the classical information in the most optimal way as the periodic 
orbits have to be calculated only once for each e. 

The symmetry reduced zeta functions corresponding to the sub-spectra with 
fixed parity tt z = ±1 are obtained by restricting the classical dynamics to the 
fundamental domain, see Appendix A. The symmetric and anti-symmetric zeta 
functions differ in phase factors, which are due to the different boundary conditions 
as explained in detail in Appendix A. These phase factors are included here in 
the Maslov index a, which can be read off directly from the three letter symbolic 
dynamics introduced in section 1.1. We obtain [31] 

er p = 2n p + n 2 + (1 - 7r z )ni, (34) 
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where n p denotes the length of the symbol string corresponding to the periodic 
orbit, and n\, n 2 count the number of symbols T, '2' in the code. Note, that the 
periodic orbit 2 running along the symmetry line has to be excluded in (33) for 

7T Z = -1. 

The product representation of the zeta function (33) is absolute convergent only 
for Im(s) > h a > [47]. For systems with a well defined symbolic dynamics 
a so called cycle expansion of products like (33) has been proposed [28, 29] to 
provide a representation with larger analyticity domain. The product is expanded 
by multiplying out the single factors and regrouping the terms in such a way that 
maximal cancellations occur. The ordering scheme makes intensive use of the self- 
similar structure of the dynamics reproduced by the symbolic description of the 
flow and groups together orbits and pseudo-orbits into contributions of increasing 
total symbol length, the so called curvature terms. For a symbolic dynamics with 
a three letter alphabet {0,1,2}, this has the form [29] 

oo 

r 1 = n( 1 -^) = 1 -E c " ( 35 ) 

V n—1 

with 

Ci = [to + ti + £2] (fundamental term) 

C 2 = [(£01 - t tl) + (t 02 ~ t t 2 ) + (t 12 - M 2 )] 

c 3 = [(^001 — ^0l) + (^011 — ^Ol^l) + (^002 — ^02) 

+ (^022 — ^02^2) + (^112 — ^12) + (^122 — ^12^2) 

+ (^021 — ^02^l) + (^012 — ^0^12 — ^01^2 + ^O^l^)] 
C 4 = . . . . 

Note that each term in the curvature contributions c n is accompanied by a 'shad- 
owing' term differing in sign. The only exception is the leading term c\, which is 
also called the fundamental term [29] in the expansion. 

The method has been applied successfully for variety of strictly hyperbolic sys- 
tems [48, 29, 49, 50] and yields an analytic continuation of the zeta function in a 
strip h a > Im(s) > h c , h c < with exponentially decreasing curvature terms c n in 
this regime. Analytic semiclassical expressions below the critical value h c can be 
obtained by modifications of the Gutzwiller-Voros zeta function itself [51], which 
will not be discussed here. Hyperbolic systems are characterized by Lyapunov expo- 
nents X p = of periodic orbits strictly bounded away from 0. The DKP, however, 
is an intermittent system even for scaled energies e above the critical value e c and 
we obtain A ~ — >■ for periodic orbits with an increasing '0'- tail in the sym- 
bol code. This behavior has strong implications on the curvature expansion and a 
possible analytic continuation of the zeta function, which will be discussed in the 
next paragraph. 

3.1 Cycle expansion with infinite symbolic dynamics 

The dynamics of the DKP is strongly influenced by the regular behavior in the 
asymptotically separable channel as discussed in detail in section 1.2. Periodic orbits 
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with equal symbol length contribute very differently to the zeta function depending 
on the time spend far from the nucleus. The weight t p of channel trajectories 
drop off algebraically with the period, whereas contributions from periodic orbits 
localized in the center region decrease exponentially in amplitude with increasing 
symbol length. As a consequence, shadowing properties in a cycle expansion (35) 
between orbits with equal symbol code length are poor. Furthermore, the regular 
orbits dominate the curvature terms c n with increasing symbol length n giving rise 
to an algebraic decay in the curvature expansion 

C n « n -5/6 e 27ri Se » 

for real s. As a consequence, an expansion of the zeta function grouping together 
single orbit terms like in (35) diverges for Im(s) < 0, i.e. exactly in the region, 
where the quantum resonances are expected. 

Controlling the divergences introduced through the coexistence of regular and 
chaotic motion is one of the main challenges in classical and semiclassical periodic 
orbit theory [52]. In the following, we present a method which overcomes the prob- 
lem of including intermittency originating from one marginally stable orbit in the 
cycle expansion of a zeta function. We point out, however, that the distribution 
of marginally orbits in a generic system with mixed regular and chaotic motion, as 
e.g. in the classical DKP for e < e c , may itself have a complicated structure. 

The marginally stable behavior for \z\ — > oo enters the symbolic dynamics 
through the symbol '0', which corresponds to half an oscillation perpendicular to 
the magnetic field axis. The symbols T and '2' are always related to motion in 
the near core region. The symbol '0' can occur everywhere in the regular channel 
arbitrary far from the origin. A single '0' in the symbol code contributes differ- 
ently to the action or stability exponent of a periodic orbit depending on how far 
from the nucleus the corresponding oscillation takes place. The classical dynamics 
is thus not reflected by the symbol '0' alone, but by the number of '0's occurring 
successively in a row. 

Consequently, we have to change our 3 - symbol alphabet to a double infinite 
alphabet in the following way; substrings starting with a T' or '2' and followed by 
n — 1 symbols '0' are transformed according to the rule 

lO" -1 -^ 1 , 20" _1 ^n 2 ; n = 1, 2, . . . , oo. (36) 

The index n corresponds to the total length of the substring in the ternary alphabet. 
The symbol code length of a periodic orbit in the infinite alphabet is given by the 
number of symbols differing from '0' in the ternary code. All periodic orbits in the 
families 10" _1 , 20 n ~ 1 correspond now to orbits n 1 , n 2 with symbol length one in the 
new alphabet. The orbits with symbol length two can be written as 

(run) ,(mn) ,{mn) m, n = 1, 2, . . . oc. 

The upper index denotes here the starting symbols T' or '2' of the different '0' 
substrings, the lower index corresponds to the length of the substrings. We realize, 
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Table 1: Infinite alphabet compared with the ternary alphabet for the shortest peri- 
odic orbit families; the last columns denote the leading coefficients in the expansions 
(20), (23) for e = 0.5. 
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that the new symbolic dynamics has a matrix structure, e.g. (mn) 11 corresponds 
to all codes lO^lO"" 1 , (mn) 12 includes all strings 10 m " 1 20 n - 1 and so on. Table 1 
lists the possible symbols up to length three in the new alphabet together with some 
periodic orbit families. A periodic orbit family is defined in analogy to section 1.2 
as all periodic orbits which differ only in one number in the lower index. Note that 
a cyclic shift both in the upper and lower index denotes the same periodic orbits. 
The symbol matrices (mn) 12 and (mn) 21 are thus equivalent, i.e. (mn) 12 = (nm) 21 
and the matrices (mn) 11 , (mn) 22 are symmetric. 

The cycle expansion in our new double-infinite alphabet has now the form 



r 1 = n( 1 -^) = 1 



71 = 1 



(37) 



with 



Cl 
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c 3 



c 4 
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EE 
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1/11 .1 , 1,1,1 ,lx 
2 km b n ' Q b k b m b n) 



+ (t 



112 

kmn 

122 
kmn 



+1+12 _ 1/11 ,2 , 1-1.1 .2\ 
h k h mn r> L km L n ' c^kmn) 



+12 +2 _ 1,1,22 , 1,1,2 ,2x 
b km b n 2 k mn 2 



1/22 y.2 _|_ 1/2,2 ,2\ 
2 km L n ~ Q L k L m L n) 



oo oo oo oo 

EEEE--- 

(=1 k=l m=l n=l 



The weights t p are written here as the elements of infinite dimensional tensors 
of rank n c , where n c is the length of the symbol string in the infinite alphabet 
or equivalent the length of the upper index. The curvature terms c„ in the cycle 
expansion (37) are now sums over all the elements of the tensors of rank n including 
outer products of tensors with rank n' < n. The prefactors are introduced to avoid 
double counting of periodic orbits. Note that the sum over the prefactors alone 
(including the sign) adds up to zero in each curvature term. 

We focus first on the leading or fundamental term c\ in the cycle expansion, 
which is a sum over all members of the periodic orbit families n 1 , n 2 . In order to 
study the behavior of a single periodic orbit family sum, we introduce the notation 



i a (s) 



oo 



71 = 1 



oo 

e 2Tri(sS a (n)-n(N+l/2))-\ a (n)/2 

71=1 



(38) 



with an upper index a = {1,2}. An additional phase 2jrNn is added (with N 
integer), which do not affect the sum. Its meaning will become clear in what 
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follows. (Note there is another phases tt/2 in the case of a = '2', see also Eq. (34).) 
Using the asymptotic results (20), (23) from section 1.2, we write 



S a (n) = n(e + -n-^+S a +S a 2/3 n-^ + j:S?n-*; (39) 
jjlogn + l a + l a 2/3 



X a (n) = ^\ og n + l a +l a 2/3 n- 2 / 3 + £/>-\ (40) 



and S is written here in units of 2n. The family dependent coefficients So, S2/3, ■ ■ ■ , lo, 
are obtained from a polynom fit on periodic orbit data up to a symbol length n = 500 
in the ternary alphabet. The leading coefficients and for e = 0.5 are listed 
in Table 1. 

The sums t a diverge for Im(s) < due to the logarithmic behavior in the 
stability exponent in (40). This is the same kind of divergence discussed earlier, 
and it appears in all other curvature terms. An analytic continuation of the sums 
(38) for negative imaginary s values is provided by the following technique [32]; first 
we transfer the periodic orbit sum into a sum over integrals, using the identity 

S(x — n) = ~y exp(27rira;) (41) 

n=— oo r=— oo 

obtained by Poisson summation. The periodic orbit sums t a can thus be written as 

1 00 roo 

t a ( s ) = -t a l (s)+ Y / dxe 2n[rx t a (s,x + l). (42) 

The continuous variable x in t a (s,x) corresponds here to the discrete symbol length 
n in (38). We now fix the index iV in (38) by the condition 

Re(s) e + - > N > Re(s) e - -, 

i.e. we change iV exactly at the Landau thresholds (32) Re(s) = (N + l/2)/e. The 
r = integral can then be evaluated by rotating the line of integration onto the 
negative imaginary axis which corresponds to a transformation x — > —\x. The 
integrand is now exponentially decreasing for all imaginary parts of s and Re(s) / 
(N + 1/2) /e. The integral is thus convergent for Im(s) < 0, but we pick up an 
essential singularity at the Landau thresholds. By the same technique, we find the 
analytic continuation for the other r - integrals. The real axis has to be rotated 
according to x — >■ ±ix and the ±-sign refers to the sign of r. We end up with the 
compact expression 

1 r°° 

t a ( s ) = -t$( s )-i dxt a (s, 1 -ix) (43) 

2 Jo 

PCO g — 27TX 

+ \/ dx^—p^[t a {s,l + ix)-t a {s,l-ix)]. 



The sum over r for r ^ is nothing but a geometric series which gives rise to 
the second integral in (43). The transformation from the original sum to the two 
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integrals in (43) is indeed an equality for Im(s) > 0, but provides an analytic 
continuation for Im(s) < 0. The second curvature terms can be worked out in the 
same spirit. The leading contributions to the action and stability exponents are 

S a (n) = n + \ n ^) + m ( e + \ m ^) + S ® + 0{n~% + m~%) (44) 

\ a (n) = -log(nm) + 1% + 0{rri + m~i) 

and 'a' corresponds to '11', '12' or '22', here. The analytic continuation of the 
double sum c 2 in (37) leads to a double integral, the higher curvature terms give 
rise to higher dimensional integrals, respectively. 

Before presenting results for the full cycle expansion, we will study the lead- 
ing term C\ in detail. We will show, that this term is indeed fundamental in the 
sense, that it contains already all the information about the gross structure of the 
spectrum. 



3.2 The QEBK - Quantization 



The dominant contribution to the periodic orbit family sums (38) is contained in 
the r = integral in (42). The other r - terms give a contribution suppress by the 
factor e~ 27rx , as can be seen from Eq. (43). To obtain the dominant s-dependent 
behavior of the periodic orbit family sums, we approximate the integral 



jo( s ) = J°° dx e 2ni[sS"(x)-(N+l/2)x]-\"(x)/2 

by stationary phase. This leads to the phase condition 
j-[u~S a {x)-{N+ l -)x \ =0, 



(45) 



(46) 



where we set s = u + iv. The solutions of (46) are in leading order (up to corrections 
0(x~ 5 / 3 )) given by 



X (u) 



u 



1 3 / 2 



2(N + \ - ue) 



(47) 



Note that x tends to infinity at the Landau thresholds u = (N + |)/e. We end up 
with the approximation to the integral 



I a (s) 



d 2 S\ . 



-1/2 



3 27n \sS a (x )-(N+l/2)x -l/8 -A a (x )/2 



(48) 



Z >1 



y/3 e- l ° /2 u~ l/2 e 2,ri [ flXo(e+ 2 x <' 2/8)+flSS ~ (JV+1/2)xo ~ 1/8 ] _|_ 0(x~ 2/3 ) 



and x (u) is the solution of Eq. (46). Note that the exponent 5/3 obtained from the 
leading behavior of the stability exponents in (40) exactly cancels the singularity 
introduced through the second derivatives of the actions at the Landau thresholds. 
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The additional term 1/8 in the exponent is due to the stationary phase approxima- 
tion. 

In a next step we approximate the term c\ in the cycle expansion (37) by the 
dominating periodic orbit family only, i.e. by the family with the smallest coefficient 
l in (40). For e = 0.5, this is the '1' - family with l = 1.696 . . ., see Table 1. This 
leads to the approximate quantization condition, 

C 1 (s)*l-i 1 (s) = 0. (49) 

Inserting the stationary phase expression (48) in (49), we can write down a quan- 
tization condition for the real parts of s directly, 

uS l (x (u))-(N + l)x (u) = M + l (50) 

Z o 

with M integer. The quantization conditions (46) together with (50) can be inter- 
preted in terms of the semiclassical Einstein-Brillouin-Keller (EBK) quantization 
of integrable systems, (see [22] for an overview). The Hamiltonian of a classical in- 
tegrable systems with / degrees of freedom can be written in the form H = H(J), 
and Ji, i = 1,...,/ denotes the conserved momenta (or actions). An example 
of such an Hamiltonian is given in (11). The EBK treatment leads to the simple 
quantization condition Jj = rij + f and the energy eigenvalues are obtained from 
inserting the quantized actions in the classical Hamiltonian. The integer numbers 
Vi are called the Maslov indices and correspond to the number of caustics along a 
path ip(t) = u){t on the torus. 

For the integrable Hamiltonian (11), we can express the action after one oscil- 
lation along the magnetic field direction by the equation (15), i.e. 

J{a,E) = J c {a,E) + aJ h {a,E), (51) 

The variable a corresponds to the winding number in Eq. (16). Due to the condition 

Mc = _ a dJ K 

da da 



" — -a valid for fixed energy, we can write the additional condition 



T a = ^ < 52 > 

for E constant. Equations (51) and (52) are indeed equivalent to (50) and (46), if 
we make the natural identification J = uS 1 and a = x and choose the quantization 
conditions 

J c (s,e) = M + ^ ,J h (s,e) = N + ^. (53) 

The harmonic oscillator is quantized with the usual term 1/2 corresponding to an 
integer Maslov index 2. We are lead to a non integer Maslov phase 1/2 or a term 
1/8 in the quantization rule for the action J c in the Coulomb direction. (Similar 
results have been found in Helium, see Ref. [32]). A full winding in the J c coor- 
dinate implies a "visit" of the near core region, where the torus approximation of 
the full dynamics breaks down. This fact may motivate the unusual Maslov term, a 
clear understanding within the EBK - theory is, however, still missing. Equations 
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similar to (51-53) have been applied to quantize stable islands in the classical phase 
space beyond the harmonic oscillator approximation assuming near-integrability of 
dynamics in the neighborhood of the stable fix-point [53, 54]. Note, that the quan- 
tization rules are derived here directly from Gutzwiller's periodic orbit formula and 
provide a quantization both in the real and imaginary part of the s - variable. 

An effective quantization scheme in terms of equations (50) and (46) is obtained 
by first quantizing the winding number x = x^,m according to 



L as 



T dx 



1, ,„ 1 



S{x){N + -)-— {N + -)x + (M + -) =0. (54) 



2 



The eigenvalues un,m are finally given by the formula 

_ (N + 1/2) x N , M + (M + |) 

un,m ~r. (obj 

j(Xn,m) 

Note that there is no free parameter. The action S(x) is determined by Eq. (39), and 
the family dependent coefficients in the expansion are obtained from the periodic 
orbit data. 

Asymptotic solutions of the Eq. (54), (55) near the Landau thresholds can be 
given in analytic form, 



L N,M 



N + \ 2(7V + I) 



■ M + l/8 _ v 1 N + 
AT + 1/2 o{ ' 2e 2 M + l/8_ 



M> 1,(56) 



where the parameter Sq (e) denotes the leading coefficient of the expansion of the 
action Eq. (39) of the dominating periodic orbit family, (which is the n 1 family for 
e = 0.5). 

In Tables 2-4, we list the energy eigenvalues obtained from our EBK - like quan- 
tization (54) and (55), which we call quasi-EBK (QEBK) quantization from now 
on [32] . They are compared with the real parts of quantum eigenvalues calculated 
from Schrodinger's equation directly as described in section 2. Our simple quasi - 
integrable approach can reproduce the quantum spectra already surprisingly well. 
The QEBK results deviate in large parts of the spectrum not more than 10% of 
the mean level spacing and the bound states in the series below the first Landau 
threshold (Re(s) < l/2e = 1) are completely reproduced within an error of S < 
4%. This makes it possible to assign approximate quantum numbers (TV, M) to 
individual states in the spectrum. Note, that there are more quantum resonances 
than approximate quantum numbers obtained from our quasi separable ansatz. The 
eigenstates at Re(s) = 2.44, 3.73 and 4.64 in Table 3, 4 cannot be reproduced by the 
QEBK - approach. These states have thus no equivalent in the separable problem 
(11), and are an effect of the coupling region in the full potential V(p,z) in (2). 
In other words, the phase space volume of the non-separable problem is (for finite 
\z\) larger than the corresponding volume in the separable approximation, which 
causes additional resonances. Their number can be estimated by the difference in 
the phase space volume and is expected to increase like s 2 in the spectrum. 
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1.0 
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Table 2: The spectrum of hydrogen in a constant magnetic field for m? = + and 
e = 0.5; the 3 rd and 4 th columns contain the quantum results for s = calculated 
by the complex rotation method, see section 2; the approximate quantum numbers 
N and M are obtain from the QEBK results in column 5. The zeros of the cycle 
expanded zeta function are presented in columns 7 and 8; the deviation 8 of the 
real parts of the semiclassical and quantum results (in percent of the mean level 
spacing in each Rydberg series) is given in columns 6 and 9. 

The approximate quantum numbers also help to clarify the structure of the spec- 
trum. The Rydberg series converging to the different Landau thresholds Re(s) = 
(N + l/2)/e overlap for N > 1, the bound state series N = is the only one remain- 
ing unperturbed. We identify 3 perturber states in the N = 1 series and 9 perturber 
in the N = 2 series. The number of perturber states increase rapidly with increasing 
N. The perturber states interfere with neighboring Rydberg states causing level 
repulsion both in the real and imaginary part of the eigenvalues. The interference 
effects are a result of the non-separability of the quantum Hamiltonian and can not 
be resolved within the QEBK approach. As a consequence, the error in the QEBK 
eigenvalues is maximal in the neighborhood of a perturber and a unique assignment 
of quantum resonances by approximate quantum numbers becomes questionable 
already above the third Landau threshold. 

We conclude, that the gross structure of the spectrum is reproduced by our 
new QEBK - quantization. To resolve the fine structure in the spectrum and the 
additional resonances, we have to include the chaotic part of the dynamics. This 
can be done in a systematic way by performing the cycle expansion (37) including 
higher order curvature terms. We will present results in the next section. 

The QEBK- approach is in particular useful to derive analytic expressions like 
equation (56) close to the Landau thresholds. In order to demonstrate the accuracy 
of the QEBK formulas in the large M limit, we insert the real part of the quantum 
eigenvalues from Table 2 - 4 on the right hand side of (56) and replace Sq on the left 
hand side by a variable Sn,m to fulfill the equality sign. The semiclassical theory 
then predicts Mium^^ Sn,m equals a constant, which is independent of N and the 
limiting value is given by the classical coefficient Sq = —0.6483. Results for different 
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Table 3: The same as Table 2 for 1 < Re(s) < 3. 
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Table 4: The same as Table 2 for 3 < Re{s) < 5. 
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J_l 

20 

Figure 5: The quantity Sn,m for the different N - Landau channels versus the M 
- Rydberg quantum number for the m 7T = + and e = 0.5. The dashed dotted line 
corresponds to the semiclassical limit Sq = —0.6483 . . .. 

Landau channels and e = 0.5 are presented in Fig. 5 and support our conjecture in 
detail. 

We conclude that the spectrum of the DKP is purely integrable deep in the 
Landau channel with two important modifications; the Maslov phase for the mo- 
tion along the magnetic field axis is non integer and the spectrum depends on an 
additional parameter <Sg( e ) which can be determined by the classical dynamics. 

The QEBK-approximation allows for a quantization of both the real and imag- 
inary part of s. The quantization condition for Im(s) can be obtained by inserting 
the quantized winding numbers xn,m an d energies u^,m in (48) and (49) which 
yields the asymptotic expressions for v = Im(s), 

A decrease in the width of the resonances proportional to M~ 3 is indeed typical for 
atomic systems [55]. 

The antisymmetric spectrum tt z = — 1 is obtained by a sign change in (49) or 
equivalent by replacing the phase factor 1/8 by 5/8 throughout this section. The 
comparison between the QEBK quantization and exact quantum calculations yields 
results of the same quality as for the tt z = +1 spectra. 

3.3 Cycle expansion results 

So far, we can describe the gross structure of the spectrum in terms of a quasi 
integrable approximation, which can be derived from the fundamental term c\ in the 
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n = 1, . . . oo 

m < 3 
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111, 112, 122, 222 


(n, 


m, k); 


n = 1, . . . oo 
m + k < 3 
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1111, 1112, 1122, 
1212, 1222, 2222 


(n, 


m, k, I); 


n = 1, . . . oo 
m + k + / = 3 



Table 5: Cutoff of the tensor summations in the cycle expansion (37). 



cycle expansion (37). The finer details in the spectrum, as e.g. level repulsion due to 
the overlap and interference of different Rydberg series is a manifestation of the non- 
integrability of the system and can not be resolved in the QEBK-approach. The 
non-integrability enters the semiclassical formalism through the chaotic dynamics 
in the near core region and its coupling to the regular channel. The full information 
about this part of the classical dynamics is contained in the higher curvature terms 
c n in the cycle expansion (37). An analytic continuation of the multiple infinite 
sums below the real s axis leads to multi dimensional infinite integrals as discussed 
in detail in section 3.1. The one dimensional integrals (43) in the fundamental term 
ci can be calculated "exactly" by standard numerical techniques. A full treatment 
of the higher curvature terms demands a determination of periodic orbit dependent 
coefficient matrices in the expansion of the actions and stability exponents like 
in Eqs. (44), and an evaluation of higher dimensional integrals. At present, this 
exceeds the limits of our computational abilities. 

Instead, we approximate the higher order curvature contributions in (37) by 
carrying out the analytic continuation only in one of the summation indices and sum 
over the first few numbers in the other indices directly. (Summing over the lower 
indices m, k or / in Table 5 up to infinity leads of course to the same divergences 
for Im(s) < as described for the single sums in chapter 3.1. This is, however, 
no problem as long as m, k, I is small.) The treatment is equivalent to a cycle 
expansion in terms of periodic orbit families of the form G0 n_1 , where the head 
string G is here a substring of length no in our old ternary alphabet starting and 
ending with a T' or '2'. The results presented below belong to a cycle expansion 
including all periodic orbit families up to a maximal head length n G = 4. The 
corresponding symbols in the infinite alphabet are listed in Table 5, the leading 
order family dependent coefficients in the expansion of the actions and stabilities 
are included in Table 1. Our expansion of the zeta function thus contains 54 periodic 
orbit families and it includes contributions up to the 4 th curvature term in the cycle 
expansion (37). 

Fig. 6 shows the absolute value of the cycle expanded zeta function for Im(s) = 
together with the real parts of the quantum resonances marked as ticks on the 
Re(s) - axis. The minima of the zeta function coincide very well with the quantum 
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Figure 6: The absolute value of the cycle expanded zeta function for the rri* = + 
subspace and e = 0.5; the ticks on the s axis denotes the position of the quantum 
resonances and the arrows indicate perturber states due to the overlap of different 
Rydberg series. (The quantum numbers (N,M) are obtained from Table 2-4.) 

spectrum and the Rydberg series structure is reproduced leading to an essential 
singularity of the zeta function at each Landau threshold Re(s) = (N + l/2)/e. 
The perturber states appear as strong modulations in the Rydberg - oscillation, 
(indicated by the arrows in Fig. 6). To obtain both the real and imaginary part 
of the quantum eigenvalues, we calculated the zeros of the zeta function in the 
complex s - plane. Results are shown in Fig. 7 in comparison with the complex 
eigenvalues obtained from our quantum calculations. The numerical values are 
listed in Tables 2-4. Our semiclassical zeta function indeed contains the information 
about both the position and width of the resonances. The strong oscillations in the 
imaginary parts of the eigenvalues are well reproduced and the error is of the same 
order of magnitude as for the real part. We also resolve the first states missing in 
the QEBK approach at Re(s) = 2.443 and 3.733. Also the ground states of the 
N = 2,3,4 series perturbing the lower Rydberg series are now obtained including 
level repulsion effects! The error S in units of the mean level spacing is uniform in 
the neighborhood of these perturbing states. In addition, the overall deviation of 
the semiclassical eigenvalues from the quantum results decreases compared to the 
QEBK-approach. A remarkable exception is the series of bound states N = 0. The 
real part of the zeros of the zeta function reproduces the quantum data again well, 
but we obtain (in contrast to all the other Landau series) here a positive imaginary 
part! A correction of the systematic deviation of these zeros from the real s - axis 
are still an open question. The occurrence of bound states itself is a pure wave 
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Im s 




Re s 



Figure 7: The zeros of the cycle expanded zeta function (+) compared with the 
quantum results (o) for m 7T = + and e = 0.5 plotted in the complex s plane. 



effect, and is associated with the largest wave length possible in the system. We 
thus expect higher ft-terms to play a crucial role in the semiclassical expansion of 
the low lying part of the spectrum, here the N = bound states. Similar effects 
have been observed in scattering billiards [56, 57], where the next to leading h 
corrections have been calculated explicitly. 

Our cycle expansion of the zeta function can so far not resolve all quantum states, 
zeros are missing near strong perturbations of the Rydberg series, see Tables 3, 
4. (Note, that these perturber can still be observed as modulations in the zeta 
function, Fig. 6). The reason for not reproducing all the quantum eigenstates is well 
understood and not a principle problem of the theory. A semiclassical quantization 
of a generic system with more than one degree of freedom demands periodic orbits 
of increasing length or equivalent higher curvature terms to resolve the increasing 
level density in the spectrum or in our example the increasing number of perturbing 
resonances. The number of orbits itself grows exponentially with the period, which 
usually sets numerical limits to the maximal curvature contributions attainable [58]. 

In addition, we approximate here the curvature terms itself by the shortest 
periodic orbit families only. We miss thereby information about the dynamics in 
the classical Landau channels by periodic orbits visiting the far regular regions two 
or more times. These orbits build up perturbations of the Rydberg series near the 
Landau thresholds due to the overlap of different N - series. As a consequence, we 
can not resolve perturber states with large imaginary part in the Rydberg series 
for large M quantum numbers. Both effects can be read off from Fig. 8, where we 
plotted the modulus of the curvature contributions for Im(s) = 0. The curvature 
contributions drop off exponentially for fixed s, and we obtain a decrease over 3 
decades for Re(s) close to zero! The total fall off becomes smaller inside each N 
- interval, (but is still exponential) and decreases also globally with increasing N. 
The non uniform behavior of the various curvature contributions in each Rydberg 
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Figure 8: The absolute value of the first 4 curvature contributions c n of the zeta 
function for real s - values. 

series is due to the finite approximation of the curvature terms itself, the overall 
increase is a manifestation of the cutoff in the curvature expansion. 

We conclude, that the cycle expanded zeta function with a double infinite sym- 
bolic dynamics includes all the important features of the whole spectrum. We re- 
solve both the Rydberg - like structure at the Landau levels and the perturbations 
due to the overlap of the different N - series. The deviations from the quantum re- 
sults and the missing levels can be understood in the context of the approximations 
made in the cycle expansion. A quantization of the odd parity spectrum m n = 0~ 
is obtained by changing the Maslov index according to (34) and yields results of 
the same quality as presented above. The semiclassical quantization presented here 
thus exceeds standard quantum mechanical methods in principle, which have prob- 
lems to obtain the high lying Rydberg resonances in the continuum, see [59] for an 
alternative approach. 

Our results also shed a new light on recent speculations about the origin of the 
very narrow resonances or quasi-bound states, which have been observed in Hydro- 
gen in a magnetic field [60, 18], (see e.g. the eigenvalue at Re(s) = 2.627 in Table 3.) 
Our semiclassical treatment suggests, that these extremely small imaginary parts 
are due to pure interference effects. They can not be traced back to special features 
in the classical dynamics, as e.g. to stable classical motion, but can be reproduced 
as a collective effect in a quantization using unstable periodic orbits only. This 
is in contrast to the Helium atom, where series of quasi-bound resonances could 
be assigned to a large stable island in phase space [54]. As a consequence, these 
small resonances in the DKP cannot be expected to be correlated with the classical 
relevant lines e = const. Detailed numerical studies indeed found a more or less 
'random' distribution of narrow resonances over the whole E - 7 parameter plane 



[61]. 
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4 Conclusions 



We have presented here the first semiclassical calculations of the resonance spectrum 
in the Diamagnetic Kepler problem. A proper treatment of the system can not be 
performed without overcoming fundamental problems in the semiclassical approach 
originating from the intermittency in the classical dynamics. A regularisation of 
the divergences in the semiclassical zeta function caused by the regular part of the 
classical motion is obtained by performing the cycle expansion of the zeta function 
in terms of a double infinite alphabet. The fundamental term and the curvature 
contributions are written as infinite sums over families of periodic orbits. The 
orbits in these families approach the marginal stable fixed point in the separable 
limit of the potential. An analytically continued integral representation of the 
various curvature contributions is given. The universal asymptotic behavior of the 
actions and stability exponents in the different periodic orbit families could be 
derived analytically, subleading contributions are obtained from periodic orbit data 
directly. 

An optimized version of the complex rotation technique is presented, which allow 
for a calculation of the resonance spectrum for relatively high lying Rydberg states 
at the various Landau thresholds. 

The good agreement between "exact" and semiclassical results both in the po- 
sition and in the width of the resonances demonstrates that the semiclassical zeta 
function as the leading term in an h - expansion indeed carries most information 
about the spectrum. An approximation of the zeta function itself uncovers an EBK 
like quantization of the marginal stable fixed point at infinity. This quantization 
scheme may serve as a useful tool to obtain an overall estimate of the spectrum 
and yields surprisingly accurate results especially for the ground state series and 
for the Rydberg spectrum near the Landau thresholds. Our results suggest, that 
the extremely narrow resonances in the DKP are a pure quantum interference ef- 
fect, (which is reproduced qualitatively by our semiclassical method using unstable 
periodic orbits), and are not associated with stable classical motion. 

Our results so far may serve as a basis for further projects such as the study 
of ^-corrections. Of fundamental interest is a deeper understanding of the semi- 
classical theory, when crossing the critical scaled energy e c from above. This is no 
problem for pure quantum techniques, but complicates semiclassical methods as 
presented here significantly due to bifurcations of periodic orbits and pruning in 
the symbolic description of the dynamics. 
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A Symmetry reduced symbolic dynamics 



A three letter symbolic alphabet for the 4-disk problem and the diamagnetic Ke- 
pler problem is obtained by defining the new symbols as the relative increment of 
two consecutive bounces s t s t+ i with two disks. Our choice of a ternary symbolic 
dynamics is a symbolic description which is closely related to the C^ v - symmetry 
in the 4-disk billiard and in the DKP - Hamiltonian (4) in //, v - coordinates. The 
new symbols correspond to the symmetry operations forming the C^ v group. Or- 
bits, which can be mapped onto each other by symmetry operations, will have the 
same code, here. The symbolic dynamics thus describes the motion in the desym- 
metrized fundamental domain, (which is e.g. the area \x > v > in Fig. 2), with 
hard reflection at the boundaries. The symbolic description was first introduced by 
Cvitanovic and Eckhardt [33] in the 4-disk problem. 

We denote this symbolic dynamics G = ■ • ■ g_29-\9§9\9i • • ■ with gi 6 {0, 1, 2}. 
The disk enumerating symbols S are obtained from the symbols G by the following 
algorithm: Choose one disk as the starting disk si (e.g. si = 1) and the next disk is 
given as s 2 = Si + 1 mod 4 if g x = or g x = 1 and it is s 2 = Si + 2 mod 4 if g x = 2. 
If the symbol g\ was a 1 or a 2 we continue adding 1 or 2 to obtain s 3 depending 
on the value of g 2 . If g\ was a we change the sign of the increment for the next 
symbol and choose the decrement —1 or —2 to get s 3 . Each symbol g^ = change 
the sign of the increment for the following symbols while g% = 1 preserves the sign. 

Using this rule we obtain for the periodic orbit G = 2 the strings S = 13 or 
S = 24; from G = 1, S = 1234 or S = 1432; from G = 20, S = 1324, S = 1342, 
S = 1423, or S = 1243; and from G = 10, S = 1232, S = 2343, S = 3414, or 
S = 4121. These are the four orbits drawn in Fig. 2. 

In the alphabet G we have the same symbolic description for all periodic orbits 
which can be mapped into each other with a spatial rotation or a reflection. Cycles 
which are related to each other by time reversal symmetry alone will be represented 
by two different symbol strings. 

The symmetry of a periodic orbit is obtained by identifying each symbol with a 
discrete group operator of the C^ v group and multiply the elements of the symbol 
string describing the cycle together. We identify symbol with the reflection a on 
the (i or v axis, symbol 1 with a 7r/2 - rotation or a C symmetry, and symbol 2 with 
C 2 , a 7T - rotation of the configuration space. The total symmetry of the periodic 
orbit becomes important, when looking at the periodic orbit zeta function in its 
symmetry factorized representation. 

The zeta-function for the 4-disk system factories according to the irreducible 
representations of the group C^ v in Q = C4iC4 2 CbiCb 2 Ce- The factorization of the 
single orbit contributions can be read off from Ref. [33] . The DKP in p, z coor- 
dinates, however, has only C 2 symmetry due to the invariance of the Hamiltonian 
(2) under the transformation z — > —z. The symmetric and antisymmetric repre- 
sentation of the group C 2 correspond here to the A\ and B\ representation of the 
Gi v symmetry of the regularized Hamiltonian (4). The DKP - zeta function can 
thus be written as product £ = OiCbi- The zeros of Qa 1 yields the symmetric 
eigenstates, while Qb 1 gives the antisymmetric part of the spectrum with respect to 
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parity transformation. The periodic orbits contributing to the single factors in the 
zeta function are those in the fundamental domain. 

Symmetric states belong to a quantization of the DKP in the fundamental do- 
main with Neumann boundaries along the line z = 0, antisymmetric states obey 
Dirichlet boundary conditions there. In the zeta function £b i; a periodic orbit picks 
up a phase tt/2 every time it crosses the axis z = 0. The number of crossings is 
equal to the number of l's in the code. A periodic orbit with an odd number of 
symbols T' thus changes sign in the zeta function ( Bl . In addition, the cycle 2, 
which runs along the symmetry line z = 0, does not contribute to an antisymmetric 
eigenstate. These states have a node on the symmetry axis, the orbit 2 is excluded 
in the B\ - zeta function. 

A time shift corresponds to a cyclic permutation in the symbol code <?j. Note 
however, that a time reversed symbol string corresponds not to the time reversed 
physical orbit. 

B Well ordered symbolic dynamics 

A second useful ternary symbolic dynamics is the one which has a natural order- 
ing identical to the ordering of the folds in the stable and unstable manifolds. 
This enables us to make a bisection search in the Poincare map to find peri- 
odic orbits using the method introduced in [34]. Well ordered symbolic dynamics 
W = ■ ■ ■ W-2W-iW WiW2 • • • for the 4 disk system is defined in Ref. [37]. The sym- 
bol w t G {0, 1, 2} is obtained by the increment s t +i — s t — 1 for odd time t and the 
decrement s t — s t +i — 1 for even time t. Reading the string with t > as the number 
7 = E™i w t /2 t gives the relative position of the fold in the stable manifold struc- 
ture and equivalent for negative time giving the position in the unstable manifold 
structure. 

We use this symbolic position for points on a Poincare plane to search for the 
periodic orbits. We have been able to find orbits up to symbol length 1000 with 
this method and we have for each family of orbits calculated all orbits up to length 
40 and some orbits up to length 500. 

Note, that for the well ordered symbols w t a time shift is not a simple shift 
operation in the symbol string because the symbols depends on whether the time 
is even or odd. 

C The trace of the stability matrix in the clas- 
sical Landau channel 

The stability matrix M, Eq. (21), of a periodic orbit going far out in the classical 
Landau channel \z\ ^> p in (2) can be split into two parts; by fixing an arbitrary, 
but large \zo\, we write 

M = M M sep , (58) 
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where the matrix M sep describes the linearized motion in the neighborhood of a 
periodic trajectory in the regular channel from a starting point (p, z ,p p ,p Zo > 0) 
to the end point (p' , z ,p' p ,p' ZQ ~ — p Zo < 0). The chaotic part of the dynamics 
in the near core region is contained in Mo- The matrix elements of Mo approach 
constants, when increasing the number n of oscillations in the regular channel for 
periodic orbits in the same families. The dominant contribution to the stability 
matrix M sep in the regular region can be described by the linearized dynamics of 
the separable Hamiltonian (9). Transforming our coordinate system from phase 
space coordinates (p, z) to action angle variables (J,(p), we can write M sep as 



M, 



T' 1 M, 



, T, 



'-sep ±,J -sep 

where T, T' denotes the transformation matrix 



d(<p h ,<Pc,Jh,Jc) 
d(p,z,p p ,p z ) 



(59) 



(60) 



at the start and ending points. The matrix M sep is the Jacobi matrix of the flow 
with respect to action-angle variables and has now a particular simple form, 



M sep {t) 



( 1 



V o 






1 



with a r 



d 2 H 



sep 



d.P c 



-3J; 



(61) 
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and t c is the time spend in the regular channel. Note, that J c and t c ~ Jl depends 
through (13) on the energy deposed in the Coulomb motion. The transformation 
T can be given as 



/ d(p h /dp 




V 



cop 






dtp c /dz 


J?/z* 



dip h /dpp 








d(p c /dp z 




\ 



J 



(62) 



Using the condition det T = 1, we can write down a similar expression for the inverse 
of T. The derivatives of the angle variables approach finite values for periodic orbits 
in the same family and n — > oo. The explicit expressions for the derivatives are not 
needed here. 

Inserting (61) and (62) in (58), we can approximate the trace of M as follows, 

TrM^ k 1 J*t c + k 2 Jc +k 3 . (63) 

All the expressions, which approach constants in the limit t c —> oo or E c —> 0, are 
collected in hi, k%. The time variable t c is proportional to the winding number 
a and n, and using Eqs. (16), (19), we obtain J c ~ n 1//3 . The trace can thus be 
written as 



TrM = k in 5/3 + k 2 n + 0{n 2/d ) 



2/3n 



(64) 



The coefficients ki depend on the history of the periodic orbit before entering the 
regular region and on the scaled energy e. The trace of M is independent of 7 after 
one period of the periodic orbit. 
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